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We report an upper limit on the fraction of V + A current, fv+A, in top quark decays, using 
approximately 700 pb~* of pp collisions at y/s=1.96 TeV acquired by the upgraded Collider Detector 
at Fermilab. For the decay t — > Wb — > Ivb (where f = e or /j), the invariant mass of the charged 
lepton and the bottom quark jet is sensitive to the polarization of the W boson. We determine 
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at the 95% confidence level, which represents an improvement by a factor of two on the previous 
best direct limit. 
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The decay of the top quark, the most massive funda- 
mental particle observed by experiment [lj, |2| , is particu- 
larly interesting as a direct probe of the charged current 
weak interaction at the highest energy scale presently 
available. In the standard model, the spin-^ top quark 
decays via the charged current weak interaction to a spin- 
1 W + boson and a spin-i b quark with a branching 
fraction above 99% and width T t = 1.4 GeV for a 
top mass of 175 GeV/c 2 . The lifetime of the top quark, 



Ti/Tt ~ 5 x 10~ 25 s, is an order of magnitude shorter 
than the typical strong interaction time-scale for binding 
of quarks into hadrons, h/ Aqcd ~ 3 x 10~ 24 s. There- 
fore, the top quark decays before hadronization, and the 
spin information is directly transferred to the decay prod- 
ucts. In the limit mj, — > 0, the pure V — A theory of 
the weak interaction predicts that the b quark has left- 
handed (—1/2) polarization (helicity) and the W + boson 
can only have either longitudinal (zero) or left-handed 
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(— 1) polarization. The right-handed (+1) polarization is 
forbidden. The fraction f° of W + bosons with longitu- 
dinal polarization is predicted at leading order in pertur- 
bation theory to be /° = m 2 /(2m 2 v + mf) = 0.70 @. 
The non-zero b quark mass and the higher-order QCD 
and electroweak radiative corrections modify these pre- 
dictions below the 1% level 0,0- However, the presence 
of non-standard-model couplings in the tWb vertex could 
significantly modify the polarization of the top quark de- 
cay products [5J, la 111 llfj ■ Previous results have either 
been limited by the small statistics of the top quark sam- 
ples mLULlELIi or have only set indirect limits 

In this Letter, we search for a V+A current in top quark 
decay, while assuming that the it production mechanism 
is in agreement with the standard model prediction. Wc 
further assume the absence of couplings from magnetic 
moment interactions in the tWb interaction, so that / 
is unchanged from 0.70 0. Then, the V + A fraction 
fv+A is related to the fraction /+ of right-handed W + 
bosons by fv+A = — and the V — A frac- 

tion fv-A = 1 — fv+A is related to the fraction / _ of 
left-handed W+ bosons by f V -A = f~/(l - /°). The 
W + boson polarization can be inferred from the angu- 
lar distribution of the charged lepton ^| in the decay 
W+ -► i+u, 

where the angle 9 is the polar angle of the charged lepton 
in the rest frame of the W + boson. The z-axis is defined 
to be the direction of motion of the W + boson in the rest 
frame of the top quark. We use the observable M? b , the 
square of the invariant mass of the charged lepton and 
the jet from the b quark, which is related to cos 8 by 

M? b ^^(m 2 -m 2 w )(l + cos9). 

The relation is exact in the limit nit, — >0. 

This search is based on a data set with an integrated lu- 
minosity of approximately 700 pb _1 acquired by the Col- 
lider Detector at Fermilab (CDF II) from pp collisions 
at y/s = 1.96 TeV. A 96-layer outer drift chamber [18j re- 
constructs the trajectories of charged particles and mea- 
sures their momenta in the region Id < 1 0] . An eight- 
layer silicon microstrip detector [20j provides precision 
tracking information in the region \rj\ < 2 to identify 
displaced vertices associated with b hadron decays. The 
entire tracking volume is located inside a 1.4 T magnetic 
field. Electromagnetic and hadronic calorimeters mea- 
sure the energies of particle showers. Drift chambers and 
scintillation counters provide muon identification outside 
the calorimeters. Gas Cherenkov counters |2lj | determine 
the luminosity. The data are collected with an inclu- 
sive lepton trigger that requires an electron (muon) with 
\r}\ < 1 and E T > 18 GeV (P T > 18 GeV/c) 0. 



We study three independent data samples enriched in 
it events. Two of the data samples are in the lepton+jets 
channel, with it— > W + bW~b events where one of the W 
bosons decays hadronically and the other leptonically. 
The lepton+jets event selection requires one isolated lep- 
ton with E T > 20 GeV, $ T > 20 GeV 0, at least three 
jets with Et > 15 GeV, and one or two ^-tagged jets. 
More details on the selection, the ^-tagging procedure, 
and the sample composition can be found in Ref. [22j |. 
We model the hard it process with the Monte Carlo (MC) 
event generator ALPGEN [23| with CTEQ5L parton 
densities and pythia [25| for hadronization, under the 
assumption that the top quark mass is 175 GeV/ c 2 . We 
simulate the detector response using GEANT |26j, |27j . For 
it production with V — A top quark decay, we estimate 
a selection efficiency, including the branching fraction, of 
Av-A = 3.4%(1.2%) for events with one (two) 6-tagged 
jets. Due to the lower average pr of the charged lepton 
for V — A, this is a factor 0.92 below the efficiency for 
V+A. 

For the lepton+jets sample with a single 6-tagged jet, 
the 6-tagged jet is from the same top quark decay as 
the charged lepton in approximately half of the it events. 
The background M 2 b distribution is a combination of 85% 
W+jets, modeled by ALPGEN Wbb, and 15% multi-jet 
events, modeled by non-isolated lepton+jets data events. 
Background-dominated data samples with only one jet 
or only two jets are consistent, in terms of both the rate 
and the shape of the M 2 b distribution, with our model of 
the background. In 695 pb _1 , we observe 304 candidates 
with a total expected background of 88+11 events. 

For the lepton+jets sample with two 6-tagged jets, the 
two possible M 2 b values of the charged lepton with cither 
the highest or the second highest Et 6-tagged jet are used 
to construct a 2-D distribution. In this way, we keep both 
the correct and incorrect combinations, and account for 
their correlation. The background is modeled by ALPGEN 
Wbb; here the multi-jet background is negligible. Non- 
uniform binning was applied in the 2-D M 2 b distributions 
in order to ensure sufficient MC events in each bin. In 
695 pb _1 , we find 75 candidates with a total expected 
background of 9+2 events. 

The third sample is in the dilepton channel, with ii^ 
W + bW~b events where both W bosons decay leptoni- 
cally. The dilepton event selection requires two identified 
leptons with opposite electric charge and Et > 20 GeV, 
$ T > 25 GeV, and at least two jets with Et > 15 GeV. 
More details on the selection and the sample composition 
can be found in Ref. |2^|. For it production with V — A 
top quark decay, modeled by ALPGEN as described above, 
we estimate a selection efficiency, including the branch- 
ing fraction, of Av-A — 0.72%, a factor 0.88 below the 
efficiency for V + A. The two possible M 2 b values for a 
charged lepton with either the highest or the second high- 
est Et jet, assumed to be produced by the fragmentation 
of the b quarks, are used to construct a 2-D distribution. 
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As we can reconstruct M? b from the top quark decay and 
from the anti-top quark decay, we make one entry for 
each charged lepton. The effect of the correlation be- 
tween the spins of the top quark and the anti-top quark 
is negligible here. Again, non-uniform binning in the 2-D 
Mf b distributions is applied. The background Mf b dis- 
tribution is the combination of three background types: 
approximately 50% from Z/j* — > with associated 

jets, 30% from W — > Iv with associated jets where a jet is 
misidentified as a lepton, and 20% from massive diboson 
pairs, WW/WZ. The Z/j* and diboson background Mf b 
distributions are modeled by ALPGEN. The misidentified 
lepton background is based on inclusive lepton trigger 
data where the second lepton is instead a jet (charged 
particle track) weighted by a probability for misidentifi- 
cation as an electron (muon). A background-dominated 
data sample with only one jet is consistent, in terms of 
both the rate and the shape of the Mf b distribution, with 
our model of the background. In 750 pb -1 , we observe 
64 candidates (12 ee, 24 fj,fj,, and 28 e/i) with a total 
estimated background of 20±4 events. 

The fraction fv+A is estimated by comparing the Mf b 
distribution in data with parent Mf b distributions for ti 
production with V — A top quark decay (fv+A = 0.0), 
it production with V+A top quark decay (fv+A = 1-0), 
and backgrounds. A binned log likelihood fit procedure 
is used to extract the parameter of interest, fv+A- We 
represent the imperfectly known accepted background 
cross section for each sample, <7b g , and the it cross sec- 
tion [29I l3(ij|. a t t, by nuisance parameters. The analytic 
expression for the likelihood for each sample, 

N 

£ = [YiP^,^) xG(a bg ,5^ bg )xG(a t t,S <Ttf ), (1) 

i=0 

is the product over all N bins in Mf b of the Poisson prob- 
abilities of observing ni entries in a given bin i, where the 
average expected bin content is /Zj, and the Gaussian con- 
straints on the estimated background and the predicted 
it production cross sections, as shown in Table [I] The fXi 
are given by: 



TABLE I: The input values for the nuisance parameters, and 
the values from the best fit to the combined samples. 



Mi = N 



(lain 



XV+A 



XV-A 



Xbg 



xv+aT v+a + x V - A T v _ A + Xb g T bl 

fv+AAy+AVtt 

&bg + Viii-Av+Afv+A + Av-a(^ - fv+A)] ' 

(1 - f V+A) Ay-AG tt 

&bg + <Ttt[Av+Afv+A + Ay-A^ ~ fv+A)] ' 

^bg 

°bg + CT t t[Av+Afv+A + Av~a(1 ~ fv+A)] 



(2) 
(3) 

(4) 
(5) 



Here, N data is the total number of observed events for the 
sample. The xv+A, xv-Ai and Xbg are the fractions of it 
production with V+A top quark decay, ti production with 
V—A top quark decay, and background, respectively. The 



Nuisance parameter 


Input (pb) 


Fit (pb) 


0« 

<7bg lepton+jets 1 6-tag 
<7bg lepton+jets 2 6-tag 
<7b g dilepton 


6.7 ± 1.0 
0.156 ± 0.017 
0.013 ± 0.002 
0.026 ± 0.006 


7.3 ± 0.9 
0.154 ± 0.016 
0.013 ± 0.002 
0.022 ± 0.006 



T v+A , T V _ A , and T bg are the probabilities for an event 
to occupy bin i of the corresponding Mf b distribution. 
Note that J^T* = 1.0. The combined likelihood is the 
product of the likelihoods of the three samples, but with 
one common Gaussian constraint on the ti cross section. 

The robustness of the fitting procedure has been tested 
with pseudo-experiments. For a given pseudo-experiment 
and a particular sample, the number of observed data 
events, N data , is distributed in three categories (it pro- 
duction with V + A top quark decay, ti production with 
V—A top quark decay, and background) as multinomial 
deviates according to their expected fractions. These 
fractions are first varied according to Gaussian distribu- 
tions to incorporate the uncertainties in the background 
and the ti cross sections. The events are generated from 
the relevant Mf b parent distribution for each category. 
The hypotheses that fv+A = 0.0, 0.1,.. .,1.0 are studied 
for 2000 pseudo-experiments for all samples combined, 
as well for the three samples separately. In all cases, the 
fit is unbiased and stable. An expected statistical uncer- 
tainty of 0.22 on fv+A is found for the combined case, 
while for the separate samples we found 0.36, 0.41 and 
0.49 for lepton+jets single 6-tagged events, lepton+jets 
double 6-tagged events, and dilepton events, respectively. 

The estimates of the systematic uncertainties on the 
measured value for fv+A are shown in Table [n] for all 
samples combined. The leading sources of systematic un- 
certainty arise from uncertainties on the measured jet en- 
ergy [3jJ, the background shape and normalization, and 
limited MC statistics. We determine all systematic un- 
certainties by performing pseudo-experiments in which 
the systematic parameter in question is varied and the 
resulting simulated data are fit to the default parent dis- 
tributions. All shifts are evaluated at fv+A = since 
none of the systematic sources are related to the struc- 
ture of the tWb vertex. 

The maximum likelihood fit for the lepton+jets sam- 
ple (single and double 6-tagged) yields a value of fv+A = 
0.21 + 0.28, including the effects of the uncertainty on the 
background and ti cross sections. For the dilepton sam- 
ple, we obtain fv+A = —0.64 ± 0.37. The probability to 
obtain a value smaller than the dilepton result is 10% for 
the hypothesis fv+A — 0. The lepton+jets and dilepton 
results are compatible at about 2.3 standard deviations. 
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TABLE II: The systematic uncertainties on the measurement 
of fv+A for all samples combined. 



Source 


Uncertainty 


Jet energy 


0.10 


Background modeling 


0.04 


MC statistics 


0.04 


Initial/Final state QCD radiation 


0.02 


Multiple pp interactions 


0.02 


6- tag efficiency (Et) 


0.02 


MC generator 


0.01 


Parton densities 


0.01 


Total 


0.12 




Mf b (GeV 2 /c 4 ) 

FIG. 2: The M$ b distribution for the lepton+jets sample with 
a single 6-tagged jet. 




"V+A 

FIG. 1: Likelihood distribution (see Equation for the lep- 
ton+jets and dilepton data samples separately and combined. 

Combining all samples, and the total systematic uncer- 
tainty from Table [nj the result for the fraction of V+A 
current in top quark decay is 

fv+A = -0.06 ± 0.25 (stat.+syst.). 

This value is in agreement with the standard model. Ta- 
bleQJsummarizes the fitted values for the nuisance param- 
eters. The likelihood distribution is shown in Fig. ^ The 
good agreement in the Mj h distribution between data and 
the best fit result for fv+A from combining all samples is 
shown in Figs. [21 and where the highest bins also 
contain overflow entries. For comparison, fv+A = 1-00 is 
also shown. 

In the absence of a signal, we evaluate an upper limit 
on fv+A using a Bayesian approach. The profile likeli- 
hood function is first determined as a function of fv+A, 
multiplied by a prior flat between 0.0 and 1.0, and nor- 
malized to yield the posterior distribution for fv+A- The 
upper limit at 95% confidence level (CL) is formed by 
integrating the posterior from zero to the value of fv+A 
that yields 0.95 for the integral. We verified that this ap- 
proach yields proper frequentist coverage for fv+A < 0.3; 



a small correction, derived by the Ncyman construction, 
would be applied to any upper limit greater than 0.3 to 
restore coverage in the region fv+A > 0.3. For the stan- 
dard model, where fv+A = 0, the median upper limit 
expected from combining all three samples is 0.38, with 
the 68% interval from 0.2 to 0.6. Combining all samples, 
we set an upper limit on the fraction of U±A current in 
top quark decay of 

f v+ A < 0.29 at 95% CL. 

This is an improvement by a factor of two on the previ- 
ous best direct limit [TJ- We estimate a ±0.07 (±0.09) 
shift in the quoted upper limit (measurement) for fv+A 
if the top quark mass is ±2.5 GeV/c 2 different from 
175 GeV/c 2 . Converting the results for fv+A to the 
fraction of right-handed W + bosons, we obtain /+ = 
-0.02 ± 0.07 and /+ < 0.09 at 95% CL. 
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FIG. 3: The invariant mass squared of the charged lepton 
and the highest Et 6-tagged jet, Mf bl , in five regions of the 
invariant mass squared of the charged lepton and the second 
highest Et 6-tagged jet, Mf h2 , for the lepton+jets sample with 
two 6-tagged jets. 
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FIG. 4: The invariant mass squared of a charged lepton and 
the highest Et jet, Mf bl , in four regions of the invariant mass 
squared of the charged lepton and the second highest Et jet, 
Mi b2 , for the dilepton sample. 
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